Modeling of pulsed laser guide stars for the Thirty Meter Telescope project 
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optics system and associated laser guide star (LGS) facility to correct for the 
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of the same average power; both the cw and pulsed lasers have the potential 
to satisfy the TMT design requirements for photon return flux. 
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1. Introduction and motivation 

The Giant Segmented Mirror Telescope (GSMT) initia- 
tive has been ranked among the highest priorities for 
US astronomy in the last two decadal surveys of astron- 
omy and astrophysics ^1 [5] • This ranking is motivated 
by the potential scientific impact of the sensitivity and 
diffraction-limited resolution that it is possible to achieve 
with a large-aperture telescope. The Thirty Meter Tele- 
scope (TMT) project has designed and is planning the 
construction of a diffraction-limited, large-aperture seg- 
mented telescope on the northern slope of the Mauna Kea 
volcano in Hawaii (characteristics of the site are sum- 
marized in [3]). For such a telescope, adaptive optics 
(AO) and an associated laser guide star (LGS) system 
are needed to achieve diffraction-limited capacity. 

Adaptive-optics systems are of increasing importance 
for current 0HH] and future [9"rrTT] ground-based astro- 
nomical telescopes. Large ground-based telescopes re- 
quire AO to correct for the image distortion induced by 
atmospheric turbulence. The use of AO requires a ref- 
erence light source, such as a bright natural star, in or- 
der to measure the effect of the atmosphere on the light 
wavefront. Often there is no bright natural star available 
within the isoplanatic angle from the astronomical tar- 
get, i.e., the angle out to which phase variations in the 
wavefront are less than one radian. In this case, a laser 
can be used to generate an artificial guide star. This can 
be done by exploiting Rayleigh-backscattered light from 
a pulsed laser [IS], or in the case considered here, reso- 
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nance fluorescence from sodium atoms in the mesosphere 
[TBI [T4] at 80 km to 125 km altitude, excited by either 
continuous-wave (cw) or pulsed lasers resonant with the 
sodium D2 transition at 589 nm. Such techniques were 
first proposed for astronomy in the 1980s [13 [H]; almost 
30 years later, laser guide star systems are available in 
all 6-10-m-class telescopes. 

Many of the key scientific questions to be addressed 
with TMT [17] — the nature and composition of the uni- 
verse, the formation of first stars and proto-galaxies, the 
evolution of galaxies and the interplanetary medium, the 
relation between black holes and their host galaxies, the 
formation of stars and planets, the nature of extra-solar 
planets, and the detection of extra-solar planets within 
the habitable zone of their host stars — require diffraction- 
limited imaging and thus a LGS-AO system. An example 
of the need for diffraction-limited imaging is observations 
of the crowded star field close to our galactic center. Such 
observations are technically challenging and astronomers 
have devised observation strategies to resolve the stel- 
lar dynamics and accurately locate and weigh the super- 
massive black hole in the center of our galaxy. Accu- 
rate determinations of the stellar dynamics have evolved 
significantly since the first survey of this star field by 
means of speckle imaging in the late 1990s [18], followed 
by Natural Guide Star-AO in the early 2000s [T5]; the 
best results have been achieved with the use of LGS- 
AO in the mid- to late 2000s [20 . Yet, there remains the 
need to improve these observations to resolve stars at an- 
gular distances closer to the galactic center. This will al- 
low testing the predictions of general relativity under the 
gravitational field of the super-massive black hole hosted 
in the center of our galaxy. This will be made possible 
by a large aperture, LGS-AO-assisted telescope, such as 



TMT. 

The TMT LGS Facility (LGSF) will initially include 
six lasers to form an asterism consisting of one LGS 
on-axis and five LGS equally spaced on a circle of 
35 arcsec angular diameter, with the laser-launch tele- 
scope located behind the secondary mirror [21] . Us- 
ing the Multi- Conjugate Adaptive Optics (MCAO) ap- 
proach, the LGSF will be able to provide diffraction- 
limited imaging in a field of view from tens of arcseconds 
to about 1-2 arcmin [22] . The LGSF design supports ex- 
pansion to up to nine lasers to form asterisms for other 
AO modes and radii from 5 arcsec up to 510 arcsec [23j . 

The original design for the TMT LGSF called for cw 
lasers (either solid-state Nd:YAG sum-frequency lasers 
or Raman fiber lasers). Subsequently, TMT has also in- 
cluded the option of using pulsed lasers based on diode- 
pumped solid-state sum-frequency Nd:YAG technology. 
Such lasers have been developed by scientists at the Tech- 
nical Institute of Physics and Chemistry of the Chinese 
Academy of Sciences (TIPC). In this article, we investi- 
gate the effectiveness of pulsed lasers for LGS, focusing 
on the TIPC laser format. We also compare the photon 
return from the pulsed laser to that of a cw laser of the 
same average power. 

The effectiveness of an LGS system in correcting the 
turbulence-induced aberrations in the image wavefront 
depends on having sufficient photon return from the laser 
beacons. The main requirements set for the first-light 
TMT AO system are the following: (1) a root-mean- 
squared (rms) wavefront error at zenith smaller than 
187 nm on-axis and 191 nm over a 17 arcsec field-of-view, 
(2) 50% sky coverage at the galactic pole, (3) operation 
up to a zenith angle of 65° , and (4) operation under at- 
mospheric turbulence conditions characterized by an at- 
mospheric length scale (Fried length tq [24 ) as low as 
0.1m in the direction of observation. The TMT wave- 
front reconstruction algorithm [25] has been tested under 
various scenarios and is expected to meet the maximum 
wavefront error specifications (see for instance [21]), with 
a significant margin for error, when using a LGS that pro- 
vides 900 photo-detected electrons (PDE) per wavefront 
sensor sub-aperture in an integration time of 1.25 ms. 
The specification of 900 PDE applies when the laser is 
pointed in the zenith direction, and includes a roughly 
factor of two safety margin to account for various effects 
that have the potential to increase the rms wavefront er- 
ror, including the effect of directing the laser at larger 
zenith angles. The photon return scales due to purely 
geometrical considerations as the cosine of the zenith an- 
gle (i.e., the inverse of the airmass); this scaling is used 
as a reference for comparison with the actual calculated 
directional dependence of photon return reported here. 
Some additional factors that affect the actual PDE count 
are the laser format (power, polarization, spectral char- 
acteristics, etc), various efficiencies in the propagation of 
the laser light to the launch telescope, atmospheric ab- 
sorption, the actual column density of sodium atoms in 
the mesosphere, and the relative direction of the local 



geomagnetic field with respect to the laser beam. Also 
important is the efficiency of the LGS detection system 
(i.e., the wavefront sensor). 

Optimizing a laser system to maximize the guide star 
photon return has been an area of ongoing research for 
many years. Reference [26) investigated the physical as- 
pects affecting the performance (i.e., photon flux return) 
of cw sodium lasers. The calculation included the impor- 
tant physical effects of Larmor precession of the sodium 
atoms in the geomagnetic field, radiation pressure (re- 
coil), and saturation of the optical transitions, as well as 
velocity-changing and spin-randomizing collisions. This 
was achieved using a density-matrix calculation with cou- 
pled velocity groups. A discussion was given of various 
mechanisms that are detrimental to photon return, i.e., 
Larmor precession reduces the efficiency of optical pump- 
ing, recoil depopulates the velocity classes that interact 
with the laser light, and optical saturation leads to stim- 
ulated emission, causing photons to be emitted along the 
propagation of the laser rather than back to the telescope. 

Pulsed-laser LGS systems, in the general case in which 
the pulses are not very short or very long, are more chal- 
lenging to model, due to the dynamical processes occur- 
ring on different time scales. Experiments and modeling 
efforts have been reported Refs. [2"7Tl29| . among others. 
However, none of the theoretical studies to date have in- 
cluded all of the effects discussed above for the cw case, 
which are all relevant to pulse lengths on a 100 us time 
scale. In addition, a full characterization of a particular 
pulsed laser depends on its pulse length and pulse repeti- 
tion rate and how these parameters compare to collision 
and relaxation rates and evolution time scales in the sys- 
tem, as well as the multi-mode and coherence character- 
istics of the laser. 

In this paper, we report the results of physical model- 
ing of LGS with pulsed lasers, including all of the above- 
mentioned effects. We focus on the TIPC pulsed laser 
and compute the expected photon return at the TMT 
site as a function of pulse length, pulse repetition rate, 
angular distance between the laser propagation and the 
direction of the geomagnetic field. As the first field test of 
the pulsed laser is planned for mid-2012 at the location of 
the Large Zenith Telescope (LZT) operated by the Uni- 
versity of British Columbia in Canada, we also perform 
simulations under the conditions at the LZT site. In ad- 
dition, the photon return of the pulsed laser is compared 
to that obtained with a cw laser of similar transmitted 
power. Briefly, we find that both the cw and pulsed for- 
mats appear to be able to meet the design requirements 
for photon return. The cw and pulsed formats provide 
comparable photon return, with the differences at the 
level of the uncertainty in the model. 

Results obtained from the numerical model also allow 
us to investigate various physical mechanisms occurring 
in the system that affect the photon return, often by aid- 
ing or inhibiting the generation of ground-state atomic 
polarization. These mechanisms can lead to interaction 
between seemingly unrelated parameters of the system; 



for example, we find that the optimal laser line width 
depends on the angle between the light propagation di- 
rection and the geomagnetic field direction, apparently 
due to the effect that each of these parameters has on 
the creation of atomic polarization. 

The paper is organized as follows. In Scc.[2]the system 
being modeled and various physical mechanisms that can 
influence the photon return are discussed. In Sec. [3] the 
details of the Bloch-equation model and the numerical 
solver are given. Section [4] presents the environmental 
data and the geometrical optics results that are used as 
inputs to the model. In Sec. [5] the method of integration 
over the spatial extent of the beam is described. Section 
[6] gives the results, and conclusions are given in Sec. [7] 

2. Description of the system 

When sodium atoms are subject to a strong circularly 
polarized light field resonant with the D2 transition, a 
number of physical processes occur that affect the effi- 
ciency at which light is spontaneously re-emitted back 
toward the light source. 

Left-circularly polarized light resonant with the D2a 
F = 2 — > F' transition group induces optical pumping 
that tends to transfer the atoms to the \F = 2, m = 2) 
ground state (Fig. [I]). This atomic ground-state polar- 
ization is beneficial for photon return in three ways: 

1. The cycling \F = 2, m = 2) -» \F' = 3, m = 3) 
transition is the strongest transition in the D2 
transition group, increasing the effective absorption 
cross section. 

2. Atoms excited on the cycling transition cannot 
spontaneously decay to the F = 1 ground state, 
where they would no longer interact with the light. 

3. Fluorescence from \F' = 3, m = 3) is preferentially 
directed along the light beam, leading to an en- 
hancement of photon flux observed at the location 
of the light source. 

Several mechanisms work to counteract the production 
of atomic polarization — the relative importance of these 
is related to the characteristic time scale of each. 

The geomagnetic field induces Larmor precession, 
which can act to wash out the atomic polarization. Be- 
cause circularly polarized light produces atomic polariza- 
tion parallel to the light propagation direction, precession 
will be induced if the magnetic field has a component 
transverse to the wave vector k (Fig. [2]). A magnetic field 
of 0.33 G, present at the Mauna Kea site, corresponds to 
a Larmor precession frequency of ~ 220 kHz or a preces- 
sion time of 4.3 us, meaning that this effect is important 
during the 120 us pulse produced by the TIPC laser. The 
detrimental effect of the magnetic field is strongly depen- 
dent on the direction of light propagation relative to the 
the magnetic-field direction. 

The atomic ground-state polarization is also relaxed 
by spin-randomization collisions, predominantly with O2 
molecules [30] - The time scale for these collisions varies 



with the O2 density, but can be estimated to be in the 
200-5000 us range [26] . This mechanism can thus weakly 
affect the dynamics during a 120 us pulse. Finally, transit 
of the Na atoms through the light beam (due to the mo- 
tion of the atoms or the motion of the beam itself) can 
also be considered to be a relaxation mechanism. The 
net effective atomic velocity can be estimated to be less 
than 370 m/s &g\ , which, for a typical beam size of 20 cm, 
corresponds to a millisecond time scale for transit relax- 
ation. This effect can thus generally be neglected for the 
pulse lengths in question, although we include it in the 
numerical model for completeness. 

The effective loss of atoms by spontaneous decay to 
the F = 1 ground state can be directly combatted by 
applying "repump" light resonant with the D2b F = 1 — > 
F' transition group to transfer the atoms back to the 
F — 2 ground state [13] (26] [30] . As the simulation results 
show, this technique is vital to obtaining sufficient photon 
return fluxes. Alleviation of the loss of atoms to the 
F = 1 state also has the indirect effect of reducing the 
effect of the magnetic field on the return flux, as this is 
one of the worst effects resulting from the destruction of 
atomic polarization caused by Larmor precession. 

Another process that tends to reduce the photon return 
flux for strong light fields is recoil of the atoms induced by 
interaction with the light. Each absorption/spontaneous 
emission cycle gives the atom a small momentum kick in 
the laser beam propagation direction; the cumulative ef- 
fect of these kicks is to create holes in the Doppler distri- 
bution, depleting the atomic velocity groups that are the 
most strongly resonant with the light (this is illustrated 
in Sec. 6 C ). Velocity-changing collisions of the Na atoms 



with N2 and O2 molecules act to reequilibrate the atomic 
velocity distribution, filling in the holes created by recoil. 
These collisions occur on a time scale of 15-300 us [26], 
depending on gas density, and so can act during a single 
laser pulse. 

3. Description of the model 

A. The Block equations 

In order to calculate the observed fluorescence, the evo- 
lution of the atoms is modeled using the optical Bloch 
equations for the atomic density matrix. The density 
matrix describes the magnitudes of and correlations be- 
tween basis wave functions in the ensemble of Na atoms, 
in particular, the spin wave functions of the states in- 
volved in the Na D2 transition. In order to account for 
atoms with different Doppler shifts, the density matrix is 
considered to be a function of atomic velocity along the 
laser beam propagation direction. (An additional degree 
of freedom may be included to account for the multimode 
structure of the laser line as discussed below.) The calcu- 
lation is semiclassical in the sense that while the atoms 
are treated quantum mechanically, the light fields are 
treated classically; thus, the effect of spontaneous decay 
must be included phenomenologically. Because the den- 
sity matrix describes all populations of, and coherences 
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Fig. 1. (Color online) Level diagram of the Na D2 line, showing transitions induced by circularly polarized light reso- 
nant with the D2a transition. The widths of the arrows indicate the relative transition strengths. RochesterFigl.eps 




eralization of the Schrodinger equation: 

l p = ±[H,p}+A(p)+p. 

at in 



(1) 



Fig. 2. (Color online) Precession of the atomic polariza- 
tion P, originally produced along the light wave vector 
k, about the magnetic field vector B. If the polarization 
is continuously generated and has a lifetime longer than 
the Larmor precession time, it will tend to be averaged 
about the magnetic-field direction. RochesterFig2.eps 



between, the 24 Zeeman sublevels making up the ground 
and excited states, the calculation describes, in principle, 
all saturation and mixing effects for essentially arbitrar- 
ily large optical and magnetic fields. (In practice, certain 
coherences in the system are negligible under normal ex- 
perimental conditions and can be neglected in order to 
increase the computational efficiency.) 

In order to perform numerical calculations, the veloc- 
ity dependence of the density matrix is discretized to 
describe an appropriate number n v g of velocity groups, 
each with a fixed longitudinal velocity. Because coher- 
ences between atoms with different velocities can be ne- 
glected, the complete density matrix p can be thought 
of as a collection of n vg separate but coupled density 
matrices, each of dimension 24 x 24. 

The evolution of the density matrix is given by a gen- 



Here the atomic level structure and the interaction with 
external fields are described by the total Hamiltonian 
H = Hq + He + Hb, with H the Hamiltonian for the 
unperturbed energy structure of the atom, He = — d ■ E 
the Hamiltonian for the interaction of the electric dipole 
d of the atom with the electric field E of the light, 
Hb = —fj, ■ B the Hamiltonian for the interaction of the 
magnetic moment fi of the atom with the local magnetic 
field B, and the square brackets denote the commuta- 
tor. The term A in Eq. (fTl) represents phenomenologi- 
cal terms added to account for relaxation processes not 
described by the Hamiltonian. In our case these relax- 
ation processes include spontaneous decay (omitted from 
the Hamiltonian due to the semiclassical approximation) , 
collisional spin relaxation ( "S-damping" ) proportional to 
S 2 p—S-(pS) [3T], where S is the electronic spin-angular- 
momentum operator, and the exit of atoms from the light 
beam due to motion of the atoms and the beam (transit 
relaxation). In addition, there are terms included in A 
to describe changes in atomic velocity due to collisions 
and light-induced recoil, as well as an effective relaxation 
rate for optical coherences that simulates a laser spec- 
trum with non-negligible bandwidth. These terms are 
described in more detail below. Each relaxation process 
described by A includes a corresponding "repopulation" 
process, so that the trace over the density matrix for 
all velocity groups is conserved, corresponding to conser- 
vation of the total number of atoms. The repopulation 
process describing the entrance of atoms into the beam 
is independent of p and so is written as a separate term 

Velocity-changing collisions are treated as hard colli- 
sions in which the velocity of the colliding atom is rether- 
malized in a Maxwellian distribution (no speed mem- 
ory). The internal state of the atom is assumed to be 



unchanged. 

Light-induced recoil is described phenomenologically 
by causing a fraction v r /Av v . g . of the excited-state atoms 
in each velocity group to be transferred upon decay into 
the next higher velocity group. Here v r is the recoil 
velocity and Av v . g . is the width of the particular ve- 
locity group. This model relies on the fact that v r — 
2.9461 cm/s (equivalent to a Doppler shift of 50 kHz) is 
much smaller than the typical value of Av v . g .. 

Equation pj supplies a linear system of differential 
equations for the density matrix elements, known as the 
optical Bloch equations. Thinking of p as a column vector 
of n v , g x 24 2 density matrix elements, the Bloch equations 
can be written as p = Ap + b, where A and b are a matrix 
and a vector, respectively, that are independent of p. The 
vector b corresponds to /3 and A to the rest of the right- 
hand side of Eq. (HI . 

The laser light field has a frequency component tuned 
near the D 2 a transition group, and may have an addi- 
tional 'repump' component tuned near the D 2 b transition 
group. Thus the matrix A has components that oscillate 
at each of these frequencies. Under the rotating-wave 
approximation, the overall optical frequency is removed 
from A. However, the beat frequency between the two 
light-field components remains. This beat frequency can 
also be removed from the Bloch equations in our case: 
each frequency component interacts strongly with one 
transition group and very weakly with the other, so the 
weak coupling can be neglected for each transition. If, 
in addition, the small magnetic-field-induced mixing be- 
tween the two hyperfine ground states is neglected, the 
beat frequency can be entirely removed from the evolu- 
tion equations. 

The fluorescent photon flux per solid angle emitted in 
a given direction can be found from the solution for p as 
the expectation value of a fluorescence operator [32] ■ 

B. Methods for simulation of multiple laser modes 

The spectrum of the TIPC laser is composed of three 
equally spaced modes. The proper method for including 
the effect of the modes in the model depends on the de- 
tails of the time dependence of the spectrum, including 
the coherence time between the modes and the pulse- 
to-pulse fluctuations of the mode amplitudes. Because 
these details were not known with certainty, three differ- 
ent methods were compared in the study. 

In the "coherent-modes" method, the three modes are 
assumed to be perfectly coherent and non-fluctuating. 
The atoms are subject to a field with components os- 
cillating at each of the three mode frequencies. This is 
equivalent to a single field, with a frequency equal to 
that of the central mode, amplitude modulated with a 
modulation frequency equal to the mode separation. Us- 
ing this model, the resulting return-flux signal exhibits 
beats at the modulation frequency. The width of each 
laser mode can be accounted for by including additional 
relaxation terms for the atomic optical coherences. We 
note that even when the effective line widths are broad- 



ened so that they overlap, the observed beats remain, 
because the modes are still perfectly coherent. Direct in- 
tegration of the evolution equations in this situation is 
computationally slow, because of the rapid oscillations 
in the solution. Therefore, we have implemented a Flo- 
quet technique in which the explicit appearance of the 
rapid oscillation is removed from the evolution equations 
through expansion of the density matrix in a Fourier se- 
ries. This results in a coupled system of equations for 
the relatively slowly evolving Fourier coefficients. These 
equations can be efficiently integrated to find the time 
dependence of the lowest few Fourier harmonics. 

In the "mode-hopping" method, the laser frequency is 
assumed to fluctuate between the three mode frequencies 
on a time scale much shorter than the pulse length. This 
is accomplished by writing density matrices describing 
"regions," each subject to one of the laser modes, and 
allowing atomic transit between the regions. The tran- 
sit rate corresponds to the width of each laser mode. In 
this model, the optical coherences are assumed to be de- 
stroyed upon transit between regions, meaning that the 
laser modes can be considered to be effectively incoher- 
ent. Thus no beats are observed in the return flux. 

In the "pulse-averaging" method, large pulse-to-pulse 
fluctuations in the relative mode amplitudes are assumed, 
so that during each pulse the laser power is present on 
only one of the three modes, with equal probability. This 
is accomplished by running the simulation three times, 
once with the light frequency equal to each mode fre- 
quency, and averaging the results. The modes are effec- 
tively incoherent with this method. 

C. The numerical solver 

The density-matrix evolution equations are generated us- 
ing the LGSBloch package for Mathematica, which is 
based on the Atomic Density Matrix package written by 
S. Rochester [33]. The system of ordinary differential 
equations is solved using a code written in C and based 
on the open-source ODE solver CVODE, from the SUN- 
DIALS package [34]. This solver addresses the stiffness 
of the system using the implicit backward differentiation 
formulas (BDFs) [35] ■ These methods require the solu- 
tion of a linear system of equations at each step of the 
solver; these solutions are accomplished with a precon- 
ditioned Krylov method. (In such a method an initial 
guess is improved by minimizing the residual over a sub- 
space with dimension much smaller than that of the full 
system [35] .) The rate of convergence of the method is in- 
creased by pre-multiplication with a block-diagonal pre- 
conditioner (approximate inverse of A) , obtained by set- 
ting all terms that connect density matrix elements from 
different velocity groups to zero, and then inverting the 
block for each velocity group. For the mode-hopping mul- 
timode method, the transit terms that couple the density 
matrices for the three modes are also included in the pre- 
conditioner, using an operator splitting technique that 
inverts these terms separately. For the coherent-modes 
Floquet technique, the terms for all of the included cou- 



Table 1. Standard parameters for the TIPC pulsed laser. 



Average laser power 


20 W 


Average projected laser power 




Mauna Kea 


12W 


LZT 


18.6W 


Polarization 


circular 


Number of modes 


3 


Mode spacing- 


150 MHz 


Pulse length 


120 us 


Repetition rate 


800 Hz 



Table 2. Standard parameters for Mauna Kea and LZT 
sites. The Fried length is the size of the circular aper- 
ture within which the variance of the wavefront phase 
aberrations is less than 1 radian 124 . 



pled Fourier coefficients are included in the block for each 
velocity group. 

4. Input data 

A. Standard parameters and environmental conditions 

The standard parameters used in the model for the TIPC 
laser are given in Tablefl] Environmental parameters and 
parameters specific to the Mauna Kea and LZT sites are 
listed in Table [2j The sodium density as a function of 
altitude is assumed to be described by a Gaussian func- 
tion centered at 93 km with a full width at half maxi- 
mum (FWHM) of 8 km. This is a rough estimate, as the 
measured density profile at any given time may deviate 
significantly from a Gaussian. However, our modeling 
indicates that the photon return has only a weak depen- 
dence on the width of the density profile: increasing the 
assumed width to 11km leads to a reduction in photon 
return of about 1%. For each site, the temperature and 
gas densities as a function of altitude in the sodium layer 
were obtained from the MSIS-E-90 Atmosphere Model 
[31)] ; for the Mauna Kea site, these values, along with 
the assumed Na density profile, are plotted in Fig. [3] 
Collision rates and the Doppler width as a function of 
altitude are derived from these data and the assumed 
collision cross sections given in Table [3} 

B. Laser spot image model 

To find the expected laser intensity profile, a geometri- 
cal optics simulation was performed. The laser launch 
telescope (LLT) has a diameter of 0.4 m and is nominally 
focused (absent diffraction and turbulence effects) at a 
range of 125 km. The Gaussian laser beam has a waist 
diameter (eT 2 of the peak intensity) of 0.24 m at the LLT 
pupil. 

We computed the expected laser intensity profile at 
ranges of 85, 95, 105, 115, and 125 km. This was ac- 
complished by modeling the propagation of the 589 nm 
laser light from the launch telescope for several instantia- 
tions of a single turbulence screen assumed located close 
to the LLT. The turbulence screen was modeled with 
a von-Karman power spectral density, an outer scale of 





Mauna Kea 


LZT 


Latitude 


19.8° 




49.3° 


Longitude 


155.5° 




122.6° 


Observatory altitude 


4050 m 




395 m 


B-field 


0.334 G 




0.525 G 


Na Larmor frequency 


234 kHz 




367 kHz 


B-field zenith angle 


126°18' 




160°30' 


B-field azimuthal angle 


9°45' 




16°57' 


Atmospheric transmission (zenith) 


0.84 




0.8 


Fried length (zenith) 


21cm 




5 cm 


Na column density 


4. x 


10 13 


m- 2 


Na centroid altitude 




93 kn 


i 


Na layer thickness 




8 km 





Table 3. Estimated values for cross-sections (see [26] and 
references therein; note that the factor of 1/2 described 
in [26 is incorporated into o p ^ a _ Q ). 

cr N a-N 2 7.15 x 10" 15 cm 2 
o-Na-o 2 7 x 10" 15 cm 2 
a^-o 2 2.5xl0- 15 cm 2 



turbulence [37] set to 30 m, and for various turbulence 
strengths characterized by the Fried lengths r$ = 0.10, 
0.15, 0.20, 0.25, and 0.30 m. The TIPC geometrical 
beam quality is characterized by a broadening factor of 
M 2 = 1.2. This M 2 factor is achieved in the light prop- 
agation model by introducing a slight coma aberration 
in the LLT mirror and a total wavefront error of 76.5 nm 
mis. In addition, because the LLT mirror has a fixed 
focal distance of /l — 125 km, the laser intensity beam 
profiles at altitudes lower than 125 km also include the 
effect of a focus error modeled as r 2 Ah/(2f 2 ), where r 
is the laser beam radius and Ah is the distance from the 
LGS spot to the 125 km focal distance. 

The beam profiles generated through the process above 
were interpolated as needed for intermediate values of r 
used for simulations at the Mauna Kea site. At the LZT 
site, the atmospheric turbulence is stronger, with the me- 
dian turbulence conditions better characterized by an r$ 
of 0.05 m. In this case the laser beam profile was ex- 
trapolated linearly down to ro = 0.05 m from the profiles 
available for larger values of ro- 

The LGS spot images are sampled at 0.121" and are 
64 pixels across. An example of the beam profile is given 
in Fig. H[a). Fitting the narrower dimension of the beam 
profiles with a Gaussian function yields angular sizes 
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Fig. 3. (Color online) Temperature and O2 and N2 gas 
densities as a function of altitude obtained from the 
MSIS-E-90 Atmosphere Model for the Mauna Kea site, 
along with the assumed Na density profile. (Note that 
deviations of the actual Na density profile from a Gaus- 
sian are often quite large.) RochesterFig3.eps 



ranging from 0.68" to 0.37" FWHM for r = 0.05 m to 
0.3 m. 

5. Spatial integration 

Integration of the observed return flux over the spatial 
extent (profile and path length) of the laser beam is ac- 
complished in four steps, using a method similar to that 
of [38]: 

1. Construct table of return flux vs. intensity and al- 
titude. 

2. For each pixel in the modeled beam profile, prop- 
agate light through sodium layer using the table, 
keeping track of absorption. 

3. Integrate return flux over path length. 

4. Sum over pixels. 
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Fig. 4. (Color online) (a) Interpolated beam profile from 
geometrical optics model, given in terms of the fraction of 
power in each pixel. The peak value is the Strehl ratio, 
(b) Return flux from each pixel in 10 3 photons/s/m 2 . 
Ro chesterFig4.eps 



We first generate a table of return fluxes as a function 
of light intensity and altitude. Parameters that change 
with altitude, plotted in Fig. [3j are N 2 and O2 density 
(affects collision rates), Na density (multiplicative fac- 
tor), and temperature (weakly affects collision rates and 
Doppler width). Sodium density is modeled by a Gaus- 
sian with integrated column density 4 x 10 13 m~ 2 . N2 
and O2 density and temperature are obtained from the 
MSIS-E-90 Atmosphere Model [36]. 

Both the return flux <!>(/, z) and total fluorescence into 
47r, $ tot (7, z), are obtained from the model. Return val- 
ues from the model are multiplied by the Na density and 
area of a 0.12" spot (the size of a pixel in the beam pro- 
file) at a given altitude. This results in values of fluxes 
per unit path length d£. 

Light power in the uplink beam is absorbed as it trav- 
els through the sodium layer, and reradiated as fluores- 
cence in all directions. For a sub-beam with angular size 
equal to one pixel in the beam profile, the power ab- 
sorbed per unit path length d£ is hu;d^>tot(I , z)/dt, with 



I = p ( e )/(°p£) 2 and z = £/X, where P{£) is the power at 
distance £, 9 p is the angular size of one pixel, X = sec C, 
is the airmass ((, is the zenith angle) [3S], and Uuj is the 
energy of one photon. Thus the power along the optical 
path is given by the propagation equation 

dP{£)/d£ = -hojd<5> tot (P(e)/(9 p £) 2 , £/X)/d£. (2) 

We integrate this equation for a range of initial light pow- 
ers, chosen to cover the range obtained using the given 
beam profile and average laser power. 

For the transmission of the return flux through the 
sodium layer, we assume linear absorption and so use 
the curve obtained from Eq. ^ in the low power limit, 
denoted by T\\ n {£). The observed return flux into a de- 
tector of unit area per unit path length emitted from a 
pixel-sized sub-beam at a distance £ is then given by 



d$ ohB (P ,£)/dt = J ^T Un (£)d$(P(£)/(e p £) 2 ,£/X)/d£, 

(3) 
where P(£) depends on the initial power in the sub-beam, 
Pq, and T a is the atmospheric transmission at zenith. 

Integrating d<J> obs (P > £)/d£ over the path length, we 
find the observed flux from each sub-beam as a function 
of the initial power, $ bs(-Po)- We then sum the return 
flux from all of the pixels in the beam profile. An ex- 
ample of the return flux profile is shown in Fig. R|b), 
the intermediate steps in the integration procedure are 
illustrated in Sec. 16 Dl 

6. Results 



In Sees. 6 A 6 C we illustrate various aspects of the signal 
using the return flux calculated from a particular location 
in the sodium layer, i.e., not integrated over the spatial 
extent of the beam. Sec. |6D| illustrates the steps taken 
in the spatial integration procedure, Se c. |6 E| presents the 
integrated return flux results, and Sec. |6F| presents plots 
that relate to optimization of the pulse length, repump 
power fraction, and laser line width. 

A. Time dependence 

Examples of the time dependence of photon flux return 
during a 120 us square pulse using the three different 
methods for treating the multimode laser spectrum (co- 
herent modes, mode-hopping, and pulse averaging) are 
shown in Fig.p] In Fig.pTa), no repump light is applied. 
Following the initial rapid rise in flux at the beginning of 
the pulse, the flux begins to fall off. This is primarily due 
to pumping to the F = 1 ground state, as can be seen by 
diverting some of the laser power into the repump light 
beam [Fig. IJ>tb)], which counteracts most of the fall off. 
As throughout, unless otherwise noted, we take 10% of 
the total laser power to be repump light. We also assume 
that the repump light is generated by amplitude modu- 
lation of the D2a laser field, meaning that an additional 
10% of the laser power is lost into the additional low- 
frequency sideband that is produced. Unless otherwise 
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Fig. 5. (Color online) Instantaneous return flux dur- 
ing a single pulse as a function of elapsed time, us- 
ing three different multimode treatments, both with no 
repump light and with 10% repump. Example plots 
are for Mauna Kea, zenith angle = 30°, azimuth angle 
= 190°, iavg = 27.5 W/m 2 . The oscillations present in 
the coherent-mode treatment are too rapid to be distin- 
guished on this time scale, and so appear as solid bands. 
RochesterFig5.eps 



noted, the repump light is taken to be of the same polar- 
ization as the main light field. The effect of the repump 
light is examined in more detail in Sec. |6 C| 

The rapid beats present when using the coherent 
modes treatment can be seen more clearly in Fig.[6j which 
shows the initial part of the pulse. 

B. Pulse shape 

The results of Fig. [5] are calculated assuming a perfectly 
square pulse shape. A measurement of the actual pulse 
shape produced by the TIPC laser is shown in Fig. [7] 
This shape differs from a square pulse in that there are 
finite rise and fall times, and also short spikes in the ini- 
tial part of the pulse. To investigate the effect of these 
differences, we defined an approximate analytical model 
for the measured pulse shape. The model represents the 
initial spikes as rapid exponentially decaying oscillations. 
While not a completely faithful reproduction of the mea- 
sured pulse shape, the accuracy of this representation 
should be enough to indicate whether the initial spikes 




Fig. 6. (Color online) An expanded view of the initial 
part of the coherent- modes plot from Fig. [5jb). The os- 
cillations are due to interference between the laser modes. 
RochesterFig6.eps 



have the potential to cause a significant deviation from 
the results obtained using an ideal square pulse. 

To model the measured pulse shape, we used a func- 
tional form with the pulse rise given by 

1 - e -('-'«)/^ + a osc e-^- to ^ T ^ sinVi/2), (4) 

where t is the pulse start time, r r is the rise time, a osc is 
the amplitude of the initial oscillations, r osc is the decay 
time of the oscillations, and ui is the oscillation frequency. 
The pulse fall is somewhat faster than exponential, and 
so is modeled by 



expt-Kt-tO/T/H. 



(5) 



where t\ is the time at the beginning of the fall, r/ is the 
nominal fall time, and p ~ 1.5. 

Appropriate parameters were chosen by eye to match 
the supplied measured plot, as shown in Fig. [7J Also 
shown in the figure is a pulse shape obtained by averaging 
out the initial oscillations of the formula Q. 

The calculated responses to a square pulse, a pulse 
with initial oscillations, and a pulse with the initial oscil- 
lations averaged out are shown in Fig. |8ja), (b), and (c), 
respectively. The pulses are each 120 us from the begin- 
ning of the pulse to the start of the fall, with the other 
pulse-shape parameters chosen to match those given in 
Fig. [JJ The three types of pulses are each normalized to 
have the same total pulse energy. A single-mode laser is 
assumed for these plots. 

The time-averaged flux from the more realistic pulse is 
higher than that from the square pulse, apparently be- 
cause the pulse energy is spread out over a longer effective 
time due to the pulse fall time. Averaging out the ini- 
tial oscillations has almost no effect on the time-averaged 
flux. 

These results indicate that the deviations of the mea- 
sured pulse shape from a square pulse have an insignif- 
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Fig. 7. (Color online) Laser pulse shapes: (top) The sup- 
plied measured pulse shape; (bottom) the pulse-shape 
model described in the text (solid line) and the same 
model with the initial oscillations averaged out (dashed 
line). The parameters used here are to = 22 us, r r = 
19 us, a osc = 1.9, t osc = 12 us, u — 2tt x 0.6 MHz, 
t\ = 202 us, Tf = 9.6 us, and p = 1.5. The overall ampli- 
tude is scaled to match the measurement. The nominal 
pulse length is t\ — to = 180 us, longer than the standard 
pulse length of 120 us used in the simulations. Rochester- 
Fig?. eps 



icant effect on the photon return, given the current un- 
certainty in the environmental and laser parameters. We 
therefore will use an ideal square pulse in the following 
simulations. 

C. Atomic-velocity dependence 

The return flux as a function of the atomic velocity group 
near the beginning of the laser pulse is shown in Fig. 
[9] The response to the three laser modes is seen, with 
the shape of each peak determined by the upper-state 
hyperfine structure of the D2 line. 

To study the effect of various mechanisms on the dy- 
namics of the atomic system, it is helpful to plot atomic 
populations as a function of velocity group at various 
times during the pulse. Figure [10] shows populations of 
the F — 1 ground state (bottom row), F — 2 ground 
state (middle row) and total ground-state and upper- 
state populations (top row), as a function of longitudinal 
atomic velocity using the mode-hopping treatment of the 
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Fig. 8. (Color online) Return flux in response to three 
different pulse shapes: (a) ideal square pulse, (b) pulse 
with initial oscillations and finite rise and fall times, and 
(c) pulse as in (b) but with initial oscillations averaged 
out. A single-mode laser is assumed. The nominal pulse 
lengths are 120 ps; the other pulse-shape parameters are 
as given in the caption to Fig. [7J RochesterFig8.eps 



laser modes. The left-hand column shows the popula- 
tions near the beginning of the pulse. A small amount 
of the total population has been transferred from the 
ground state to the upper-state, and some atoms have 
been transferred from the F = 2 ground state to the 
F = 1 ground state by optical pumping. The middle 
column shows the populations at the end of the pulse if 
no repump light is used. In this case the majority of the 
F = 2 atoms have been pumped into the F = 1 ground 
state, where they do not interact with the laser light. As 
a result, the number of atoms in the excited state (and 
hence the photon return flux) is greatly reduced. The 
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Fig. 9. (Color online) An example of return flux as a func- 
tion of the longitudinal atomic velocity (given in terms 
of the corresponding Doppler shift), at a time near the 
beginning of a laser pulse. RochesterFig9.eps 



hole-burning effect of atomic recoil can also be seen in 
the total ground-state population. The right-hand col- 
umn shows what happens if repump light is used. In 
this case, the population of the F = 2 state is largely 
preserved, as atoms pumped to the F — 1 state are re- 
pumped to the F = 2 state. When repump light is ap- 
plied the effect of recoil becomes much more dramatic, as 
each atom is able to interact with the light many more 
times during the pulse. 

D. Spatial integration 

Here we illustrate the steps taken to perform the inte- 
gration of the photon return over the spatial extent of 
the beam. First, we run the model a number of times 
assuming a 0.12" pixel-sized spot to construct a table of 
the return flux d$>(I, z)/d£ as a function of altitude and 
intensity (Fig. [ilk) . At the same time we also construct a 
table of the total fluorescence (Fig. [lib ) , which is used to 
calculate the resonant absorption of light in the sodium 
layer (Fig. p). 

Accounting for resonant absorption on the uplink and 
downlink, we can then construct a table of observed flu- 
orescence at the detector per unit path length in the 
sodium layer, d$> bs(Po,£)/d£, as a function of distance 
and initial power in each pixel (Fig. 13). Integrating over 
distance, we find the total return per pixel as a function 



of initial power, $ bs(-Po) (Fig. 14). As seen in Fig. 14 
the efficiency of return is constant for low power (the lin- 
ear regime), begins to rise as power is increased, due to 
the beneficial effects of optical pumping, and then falls 
off due to saturation effects, pumping to the F = 1 state, 
and recoil. 

We then take the beam profile obtained from the geo- 
metrical optics simulation with the appropriate adjusted 
Fried parameter (Fig. Ek), and apply <I?obs(-Fo) to each 
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Fig. 10. (Color online) Populations of the F = 1 ground state (bottom row), F = 2 ground state (middle row) 
and total ground-state and upper-state populations (top row), as a function of longitudinal atomic velocity for the 
mode-hopping method. Populations are shown for times near the beginning of the pulse (left column), and at the 
end of a pulse without repump light (middle column), and at the end of a pulse with repump light (right column), 
superimposed over the unperturbed equilibrium populations (dashed lines). RochesterFigl0.eps 



pixel to obtain the return flux spot profile (Fig.Hb). Due 
to the saturation effect, the central peak of the return flux 
profile is compressed relative to that of the beam profile, 
increasing the relative prominence of the dimmer features 
due to the astigmatism seen at the top of the profile. 

Summing over the pixels of the photon return flux pro- 
file gives the final calculated value for the total return 
flux. 

E. Beam-integrated return flux 

1. Return flux as a function of zenith angle 

The total return flux integrated over the spatial extent of 
the beam is shown in Fig. [l5]as a function of zenith angle 
for pulsed light, using each of the three multimode meth- 
ods, and also for narrow-band cw light. Positive zenith 
angles correspond to 190° azimuth, while negative zenith 
angles correspond to 10° azimuth. Results are shown 
both for no repump light and for the standard repump 
fraction (i.e., 10% of the light power in the repump beam, 
80% of the laser power in the D 2 a light, and 10% lost). 

The three methods for modeling the laser modes give 
results that vary by ~15%; this may be taken as an in- 
dicator of the level of uncertainty in the model due to 



incomplete knowledge of the characteristics of the laser. 
(By running calculations with finer grid spacings and 
stricter tolerances, we have found that uncertainties due 
to numerical approximations in the model are below 5%.) 
At this level of uncertainty, the narrow-band cw light 
provides approximately the same photon return as the 
pulsed format. These results indicate that both pulsed 
or cw light have the potential to satisfy the TMT design 
requirements as long as repump light is used, and will 
likely not satisfy the requirements if repump is not used. 

2. Reversed-polarization repump 

For technical reasons, it may be easier to produce repump 
light with opposite polarization to that of the main D 2 a 
light. Figure [16] shows the effect of reversing the po- 
larization of the repump light. Results are shown for 
a pulsed laser with no repump, with 10% repump, and 
with 10% repump with opposite circular polarization to 
the D 2 a light, using the pulse-averaging multimode treat- 
ment. Using opposite circular polarization for repumping 
incurs a 5-10% penalty, depending on zenith angle, rel- 
ative to repumping with the same polarization. Because 
this is a comparison of the same laser format and multi- 
mode model with itself, an observed difference at this 
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Fig. 11. (Color online) (a) Return flux d$(I,z)/d£/I 
in units of 10 9 photons/s/sr/m/(W/m 2 ) as a func- 
tion of light intensity and altitude. (b) Total flux 
d$tot(I,z)/d£/I in units of 10 9 photons/s/m/(W/m 2 ), 
used to calculate optical absorption in the sodium layer. 
RochesterFigl 1 .eps 
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Fig. 13. (Color online) Normalized observed return flux 
per d£ from a pixel-sized sub-beam, d$ o b s (P ,£)/d£/P , 
in units of photons/s/m 2 /m/W. RochesterFigl3.eps 
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Fig. 12. (Color online) Transmission through the sodium 
layer as a function of path length and initial light power 
in a pixel-sized beam. RochesterFigl2.eps 



Fig. 14. (Color online) Total observed normalized return 
flux $obs(-Fb)/-Fb from a pixel-sized beam as a function 
of initial launched light power in that pixel. Rochester- 
Figl4.eps 



level may be numerically significant. 



3. Dependence on spot size for LZT site 

To find the dependence of return flux on spot size for 
the LZT site, the beam profile was scaled proportionally 
from the nominal size 0.68" FWHM obtained from the 
geometrical optics model using the standard Fried pa- 
rameter tq = 0.05 m. When varying the beam size from 
1/2 to 3 times the standard beam size, the return flux 
varies from about 25% below to 25% above the value 
obtained with the standard size. 
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Fig. 15. (Color online) Beam-integrated return flux for Mauna Kea (as a function of zenith angle) and LZT (at zenith), 
for pulsed light using each of the three multimode methods and for narrow-band cw light. Results are given for both 
no repump light and for 10% of the light power in the repump beam. The vertical axis on the left side is in units of 
10 6 photons/s/m , while the axis on the right side of the plot for Mauna Kea shows photo-detected electrons (PDE) 
assuming a collection time of 1.25ms, collection area of 0.5m 2 , and 48% detection efficiency. The dashed gray line 
shows the design requirement for 900 PDEs at zenith, combined with the nominal dependance on the inverse of the 
airmass X expected from purely geometrical considerations. (The optical path length in the sodium layer increases 
with X, while the detector collection solid angle falls as X~ 2 .) The higher return at the LZT site is primarily due to 
the higher throughput of the laser launch system and the different direction of the local geomagnetic field compared 
to that at the Mauna Kea site. RochesterFigl5.eps 



F. Optimization 

In this section, we plot the integrated return flux for the 
Mauna Kea site as a function of zenith angle and various 
other parameters for use in optimization. Results for 
pulsed and cw lasers are given; for the pulsed format we 
use the pulse- averaging multimode treatment. 

Figure [18] shows a contour plot of return flux in units 
of 10 6 photon/s/m 2 as a function of zenith angle and 
pulse length, holding the repetition rate and the average 
laser power constant. The standard amount of repump 
light (of the same polarization as the T>2& light) is used. 
The plot runs over the range of pulse lengths that may 
be accessible to the TIPC laser. In this range, the plot 
shows that longer pulse lengths are generally better. 

Figure [19] shows a contour plot of return flux as a func- 
tion of zenith angle and repump power percentage for the 
pulsed and cw cases. For Mauna Kea, it appears that 
10% repump is optimal for both pulsed and cw light. 

Figure [20] shows a contour plot of return flux as a func- 
tion of zenith angle and the width of each laser mode for 
the pulsed and cw cases. The standard amount of re- 
pump light is used. The optimal width depends on the 
zenith angle, or more precisely, the angle between the 
light propagation direction and the magnetic field. When 



this angle is small, a broader line width is advantageous, 
apparently because the reduced effect of Larmor preces- 
sion means that less light intensity per velocity group is 
required to produce beneficial atomic polarization. When 
the angle is larger, narrower laser lines are optimal. For 
large angles, the optimum line width is reasonably con- 
sistent with the value of ~20 MHz obtained from the for- 
mula given in [5S] for the cw case. From these results 
we see that increasing the line width of the modes of the 
pulsed laser beyond the current assumed value of 15 MHz 
would result in reduced flux return for many observation 
directions. 

7. Conclusions 

Numerical simulations of photon return using both a 
pulsed laser with specifications provided by TIPC and a 
narrow-band cw laser have been performed for the Mauna 
Kea and LZT sites. The results indicate that, if repump 
light is used, both the cw and pulsed formats may be 
able to meet the design requirements for photon return. 
Within the level of uncertainty in the model, the cw laser 
format provides comparable return to the pulsed format. 
We have also investigated the dependence of the return 
flux on parameters such as the pulse length, the fraction 
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Fig. 16. (Color online) Beam-averaged return flux and 
PDE for a pulsed laser without repump, with 10% re- 
pump, and with 10% repump with opposite circular po- 
larization as the D2a light. RochesterFigl6.eps 
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Fig. 18. (Color online) Return flux in units of 10 6 
photon/s/m 2 as a function of zenith angle and pulse 
length for the Mauna Kea site. The repetition rate and 
average laser power are held constant as the pulse length 
is varied. RochesterFigl8.eps 
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Fig. 17. (Color online) Return flux as a function of spot 
size (FWHM) at the LZT site, using the pulse-averaging 
method. The dashed gray line shows the spot size ob- 
tained from the geometrical optics model using the stan- 
dard Fried parameter ro = 0.05 m. The beam profile 
was directly scaled proportionally smaller and larger to 
obtain the results for this plot. RochesterFigl7.eps 



of laser power dedicated to repumping, and the laser line 
width. For the case of the laser line width, we found that 
the optimal value has a strong dependence on the angle 
between the light propagation direction and the magnetic 
field direction. 

Future modeling work could use time-resolved mea- 
surements of the mode structure of the TIPC laser to 
obtain more accurate results. Another avenue of research 
could be to quantify the impact of changing environmen- 
tal conditions on the photon return, for example, the 
sodium-layer density profile, optical seeing conditions, 
and uncertainties in collision rates. 



Finally, experimental validation of the model will 
greatly aid future development. Published reports of 
comparison between theoretical models and measure- 
ments for Na LGS systems have been limited; there are 
some results for a cw LGS at the Starfire Optical Range 
[10"] and for a transportable cw LGS developed by ESO 
[4"T] . Significant information can be obtained even in the 
absence of absolute sodium density measurements: by 
measuring the dependence of the photon flux on experi- 
mentally variable parameters, such as zenith angle, light 
power, light polarization, and repump power fraction, the 
accuracy of various aspects of the model may be verified, 
providing information about the accuracy of the absolute 
photon return values obtained from the simulation. 
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